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ABSTRACT 

Nineteen carbon stars that show lithium enrichment in their atmospheres have been 
discovered among a sample of 674 carbon stars in the Large Magellanic Cloud. Six of 
the Li-rich carbon stars are of J-type, i.e., with strong 13 C isotopic features. No super- 
Li-rich carbon stars were found. The incidence of lithium enrichment among carbon 
stars in the LMC is much rarer than in the Galaxy, and about five times more frequent 
among J-type than among N-type carbon stars. The bolometric magnitudes of the Li- 
rich carbon stars range between —3.3 and —5.7. Existing models of Li-enrichment via 
the hot bottom burning process fail to account for all of the observed properties of 
the Li-enriched stars studied here. 

Key words: Galaxies: Large Magellanic Cloud, Stars: carbon 



1 INTRODUCTION 

Lithium is a sensitive indicator of nucleosynthesis and con- 
vective dredge-up in stars, as it is easily produced as well as 
easily destroyed in stellar interiors. 

It has been known for a number of decades (e.g. Torres- 
Peimbert & Wallerstein, 1966) that some Galactic asymp- 
totic giant branch (AGB) stars possess large surface lithium 
abundances, reaching values one or two orders of magnitude 
higher than the cosmic value. More recently, Li-rich AGB 
stars have also been identified in other Local Group galaxies, 
namely the Large and Small Magellanic Clouds (LMC and 
SMC) and M31. The great majority of Li-rich AGB stars 
discovered to date in the Magellanic Clouds have oxygen- 
rich rather than carbon-rich atmospheres (Smith et al. 1995) 
(hereafter SPLL95), although the opposite may be true in 
the Galaxy (Abia et al., 1993, Abia & Isern 1997). 

The lithium enrichment of the atmospheres of oxygen- 
rich AGB stars (i.e., stars with C/O < 1) is relatively well 
understood within the framework of the hot bottom burning 
(HBB) process during which the bottom of the convective 
envelope of the star reaches the top layers of the hydrogen- 
burning shell. Sackmann & Boothroyd (1992) showed that in 
the intermediate mass range (4-6M0) the Cameron- Fowler 
7 Be transport mechanism (Cameron & Fowler 1971) is very 
effective at producing 7 Li (via the decay of 7 Be) and can 
account for the super-Li-rich, oxygen-rich AGB stars found 
by Smith & Lambert (1990) in the Magellanic Clouds. 

Carbon-rich AGB stars (C/0> 1), (i.e., carbon stars), 
also display enhanced lithium abundances. Very large equiv- 



alent widths (W 67 07 ~ 10A) of the LilA6707 resonance 
doublet have been found in a few Galactic carbon stars: 
examples are WZ Cas and WX Cyg which are classified 
as N7 (C 2 2) and J6 (C 2 3") by Keenan (1993) and Barn- 
baum, Stone & Keenan (1996) respectively. However, most 
Galactic carbon stars have much weaker lithium lines with 
W 6 707 < 0.5A and some have 0.5A < W 6 707 < lA (Torres- 
Peimbert & Wallerstein 1966). In the Magellanic Clouds, on 
the other hand, carbon stars that display lithium-enriched 
atmospheres seem to be rarer and generally less extreme in 
terms of lithium abundance, compared with their Galactic 
counterparts or with the oxygen-rich AGB stars within the 
Clouds (e.g., SPLL95) mentioned in the previous paragraph. 
In the most comprehensive study yet of Li-rich AGB stars in 
the Magellanic Clouds by SPLL95, only four Li-rich carbon 
stars were recorded in the LMC and just two in the SMC. 

The study of Li-rich carbon stars present difficulties 
both observationally and theoretically. On the observational 
side, severe measurement difficulties arise from the blend- 
ing of the lithium lines with the complex array of lines 
in the red A 2 fI-X 2 E CN system. This heavy line blan- 
keting also means that equivalent width measurements are 
generally made relative to a pseudo-continuum rather than 
the true continuum. From the theoretical viewpoint, the 
coexistence of carbon and lithium in the atmosphere of 
an AGB star is difficult to achieve, at least in the hot- 
bottom-burning scenario, since the conditions required to 
ignite lithium production also lead to fast CNO processing 
at the base of the envelope, with the subsequent destruc- 
tion of the 12 C that was previously conveyed to the stellar 
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surface. Ventura, D'Antona & Mazzitelli (1999), however, 
showed that, for a certain range of stellar masses corre- 
sponding to a very narrow range in bolometric magnitude 
of —5.95 < Afboi < —5.75, carbon can survive long enough 
in the convective envelope to keep C/O > 1, despite the ac- 
tivation of the Cameron-Fowler mechanism. At the same 
time, the 13 C/ 12 C ratio increases and the star displays the 
spectral characteristics of a J-type carbon star which is also 
Li-rich. These models, however, cannot explain the existence 
of the fainter (by up to about one magnitude) Li-rich car- 
bon stars found by SPLL95 in the Magellanic Clouds, as 
well as by Abia et al. (1993) in the Galaxy. Abia & Is- 
ern (1997) proposed an alternative mechanism to HBB, in- 
volving deep mixing, operating in low-mass carbon stars 
(M<2M0) which can account for the fainter Li-rich and 
13 C-rich (J-type) carbon stars, found in the Galaxy and in 
the Magellanic Clouds. A more detailed description of the 
various models will be discussed later in the light of the new 
observations described in the present study. 



For Li-rich carbon stars, the derivation of reliable 
lithium abundances, the acquisition of adequate statistics 
of their occurrence, as well as an understanding of their ex- 
act evolutionary stage, are all outstanding and open issues 
which are important not only for late-stage stellar evolu- 
tion and nucleosynthesis, but also for modelling the present- 
day lithium abundance in our Galaxy and other galaxies 
(e.g., Travaglio et al. (2001), Ventura et al. (2000), Romano 
et al. (2001)). It must be emphasized that the study of 
lithium production in external galaxies such as the Magel- 
lanic Clouds is particularly important in this respect, since 
lower metallicity and a different evolutionary history would 
affect the yield of lithium on galactic scales. 



2 OBSERVATIONS 

Most of the carbon stars observed were selected from the 
Kontizas et al. (2001) catalogue of LMC carbon stars: SIM- 
BAD identifier - KDM. This catalogue was constructed by 
identifying the strong (1,0) and (0,0) C2 Swan bands (at 
4737A and 5165A respectively) in carbon stars seen in low- 
dispersion spectra on an objective-prism survey taken in the 
blue— green with the 1.2-m UK Schmidt Telescope. 

The spectroscopic observations were obtained at the 
AAT on 21 and 22 January 1998, using the 2dF multi-fibre 
spectroscopic facility, which enables up to 400 spectra to 
be recorded simultaneously, in a pair of identical fibre-fed 
spectrographs fitted with TEK 1024 CCDs (Lewis et al. 
2002) - see http://www.aao.gov.au/2df/. The target stars 
covered four 2° diameter fields. The spectra were centred 
on the Au = +2 Swan bands near 6200A and covered the 
wavelength range 5675-6785 A. Subsequent 2dF observations 
of more Magellanic Cloud carbon stars were taken over a 
slightly different wavelength range which included the strong 
(0, 1) Swan band at 5635A but excluded the Li I A6707 line. 

The 1200R gratings were used, yielding spectra with 
1.1A per pixel or an effective resolution of ~ 2.5A. Several 
exposures were taken for each field, typically 3 x 900 s or 
2 x 1200 s, together with offset sky, arc and fiat field expo- 
sures. The data were reduced with the A AO's 2dfdr data 
reduction package (Bailey et al. 2001) and the iraf package 
and yielded spectra with signal-to- noise ratio (S/N) usually 
around 30 per pixel. Full details of the observations are given 
by Cannon et al. (in preparation). 

Relative radial velocities of all the stars were obtained 
by using a cross-correlation technique in FIGARO. The high 
signal-to-noise ratio (S/N) usually attained for most stars 
and the large number of spectral features resulted in an in- 
ternal precision of about 1/20 pixel though with external 
accuracies closer to 1/10 pixel. 



Following the completion of a catalogue of 7760 carbon 
stars in the Large Magellanic Cloud by Kontizas et al. (2001) 
(see also Dapergolas et al. (1996)), an extensive programme 
of spectroscopy of carbon stars has been started (Cannon 
et al. 1999) using the 2dF multi-object spectroscopic facil- 
ity on the Anglo- Australian Telescope (AAT) (Lewis et al. 
2002). The subject of the present paper is the detection of 
the LilA6707 resonance doublet in the spectra of 19 out of 
the 674 carbon stars studied. This is the first time that such 
a large homogeneous spectroscopic sample of carbon stars 
has been studied in the context of lithium enrichment in 
any galaxy. Having homogeneous spectra for such a large 
sample of carbon stars, all at practically the same distance 
and with little reddening, provides a firm basis for making 
stronger statistical statements about the LMC carbon stars 
than can be done for (current) Galactic, or other Magellanic 
Cloud samples. The paper is organized as follows: Section 
2 describes the observations, Section 3 describes the basic 
characteristics of the spectra obtained, while Section 4 out- 
lines the method used to identify and measure the equivalent 
widths of the Li I A6707 resonance doublet. In Section 5, the 
photometric properties of the identified Li-rich carbon stars 
are investigated, while in Section 6, the results are compared 
with those of previous studies, and some theoretical impli- 
cations are discussed. 



3 SPECTRAL ANALYSIS 
3.1 General considerations 

The spectra have remarkable similarities due to the abun- 
dance of CN lines. They can be classified, however, in four 
main groups, which will be used to aid the analysis. These 
are: (1) typical N5-type carbon stars, with strong Swan 
bands and CN bands, (2) stars with weak Swan bands near 
6200 A due to higher temperature and/or lower carbon abun- 
dance, (3) stars with very weak carbon features, (4) J-stars 
with strong 13 C 12 C and 13 C 14 N bands as well as the normal 
12 C 12 C and 12 C 14 N bands. As noted earlier, identification of 
the J-type stars is important in the context of the theories 
of lithium enrichment in carbon stars. Details of the precise 
method used to identify the J stars are given by Morgan et 
al. (2003); but it is sufficient here to note that the method 
selects stars with a similar level of 13 C to the standard defi- 
nition of a J star, which is equivalent to 12 C 13 C < 15 (Abia 
& Isern 1997). Allocation of the other stars to the first three 
groups was made by simple visual inspection. Examples of 
the spectral groups can also be found in Barnbaum et al. 
(1996). (1: TV Lac; 2: BD+2°3336 and ZPsc; 3: HD 123821 
and HD 76846; 4: HD 10636 and EUAnd). Each observed 
star (with sufficient S/N) was assigned to one of these four 
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groups. The numbers of stars in the four groups were 486, 
92, 23 and 66 respectively making a total number of 674 
stars. 

The spectral range used in the following analysis was 
limited to 6690-6730A. Each spectrum was corrected for 
its radial velocity using a grid of interpolated wavelength 
points to allow this to be done to the sub-pixel accuracy of 
the velocities. It was then normalized using mean counts in 
bands at 6695-6703A and 6712-6720A. 

A mean spectrum was constructed for each one of the 
four spectral groups. The mean spectra of N and J stars 
(Groups 1 and 4) are shown in Figs la&ld respectively. 
The mean spectra of the carbon-weak stars (Groups 2 and 
3) are like that for Group 1 but with weaker features. The 
J-star combined spectrum is clearly different from the N- 
star spectrum because of the presence of 13 CN bands. The 
spectrum of the LMC carbon star BMB R-46 taken from 
SPLL95 is shown in Fig. lb; the similarity between it and the 
N-star spectrum of Fig. la is striking. The Galactic carbon 
star SYEri, which is classified by Barnbaum et al. (1996) as 
N5 (C2 5), is shown at a higher resolution in Fig. lc (Abia et 
al. 1993) and is again similar to the combined N-star spec- 
trum. However, there are differences: CalA6717 is strong 
in SYEri but is not apparent in either of the LMC stars 
and there are more features in the region 6702-6710A of the 
Galactic star. These differences are due in part to the lower 
resolution of the LMC spectra and also to the lower metal- 
licity of the LMC. Indeed, the strength of the calcium line in 
SYEri is consistent with solar abundance (Abia et al. 1993). 

3.2 Identification of Li I A6707 

The Li 1 resonance doublet at 6707.8A lies in a spectral re- 
gion which is heavily affected by many atomic and molecular 
features. These are mainly CN lines, but molecular lines of 
C2 and TiO are also present near the Li doublet. However, 
the C2 lines are of minor importance compared with the 
CN lines since at all optical depths the CN concentration 
is a factor of 7 to 100 times higher than the C2 concentra- 
tion (Abia et al. 1993). Moreover, since C/O > 1, CO for- 
mation is favoured against TiO. The presence of strong CN 
bands makes it difficult to use simple visual means to select 
stars with enhanced lithium, and even harder to measure 
the strength of the line. 

The method chosen here to identify the lithium line in a 
particular star was to subtract the mean observed spectrum 
(described above) of the spectral group appropriate for the 
star in question, scaled according to the strength of the CN 
bands in the individual spectrum. The scaling factor was de- 
termined by minimizing the difference between the spectrum 
and the scaled mean. The subtracted spectrum is mainly 
that of CN, but includes contributions from metal lines and 
other molecular bands. The subtraction process is designed 
to reveal only those features that deviate significantly from 
the scaled mean spectrum. Fig. 2 shows examples of this pro- 
cess for three stars: two typical N stars (Group 1), one with 
and the other without a lithium line, and a J star (Group 4) 
with a strong lithium line. The relative smoothness of the 
residual spectra beyond the bounds of the lithium line and 
seen throughout the entire upper panel shows that this is 
a valid procedure which takes account of the strong star- 
to-star variations in the overall strengths of the CN bands, 
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Figure 1. Spectra in the region around Li I A6707. (a) mean spec- 
trum of Group 1 (~N5) stars; (b) LMC star BMB R-46 from 
SPLL95; (c) Galactic star SYEri (type N5(C 2 5)) from Abia et 
al. (1993); (d) mean spectrum of Group 4 (J) stars. 



presumably depending on both temperature and carbon and 
nitrogen abundances. 

Residual spectra such as those shown as solid lines in 
Fig. 2 were constructed for every observed star. Lines in 
the residual spectra were then sought as follows: a second- 
order polynomial was fitted to the residual spectrum using 
the least-squares method in order to define the 'continuum' 
and all maxima and minima (potential emission and ab- 
sorption lines respectively) about this fit were identified. 
A gaussian, superimposed on the quadratic fit, was then 
fitted to each maximum and minimum allowing its central 
wavelength (A cen ), height and full width at half maximum 
(fwhm) to vary. The equivalent widths of the identified 
(both emission and absorption) lines (Wa) were then mea- 
sured. It should be emphasized that these equivalent widths 
should be used in a relative sense only, due to their strong 
dependence on the way the continuum is defined. Other au- 
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Figure 2. Normalized spectra of three carbon stars. The dashed 
line is the spectrum of a star, the dotted line is the mean spectrum 
of all stars in the same spectral group and the solid line is the 
difference between the spectrum and the scaled mean spectrum. 

(a) , the N star (Group 1) KDM 3863 which has no lithium line; 

(b) . the N star KDM 1792 which does show a lithium line; (c). 
the J star (Group 4) KDM 2832 which has a strong lithium line. 
The scaling factors applied to the mean spectra were 0.6, 1.0 and 
0.9 respectively. 



thors use a pseudo-continuum, so equivalent widths from 
different sources are not directly comparable. 

Clearly, a spectrum with poor signal-to-noise ratio will 
generate random 'lines'. It turns out that the number of 
detected lines starts to rise significantly once the count level 
falls below ~500 counts per pixel. Consequently, 53 stars 
with fewer than 500 counts per pixel were excluded from 
further consideration. 

The derived Wj are plotted against A cen in Fig. 3. What 
appears as a gap around W\ = is merely a consequence 
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6690 6700 6710 6720 6730 

Figure 3. Equivalent widths (Wj) of features detected in the 
residual spectra. Absorption features appear in the lower half 
of the figure; emission features appear in the upper half. All 
~W\ arc shown as positive values, the ordinate running posi- 
tively from zero in each half of the diagram. The symbols de- 
note the FWHM of the fitted gaussian line: FWHM<2.oA (plus), 
2.0A < FWHM < 5.0A(circle) and FWHM > 5.0A (triangle). For clar- 
ity, only one in three points is plotted for \ W\\ < 0.2 A. 



of the limit applied to the line search procedure. Absorp- 
tion features appear in the lower part of the figure; emission 
features appear in the upper half. The fairly uniform distri- 
bution of points with \ W\\ < 0.25A demonstrates the noise 
level of the identification procedure, which results from the 
S/N in the original spectra, the spectral mismatch between 
star and template and/or velocity errors. Most of the fea- 
tures in this region of the diagram have fwhm ~ 2A and are 
indeed unlikely to correspond to true absorption (or emis- 
sion) features given the spectral resolution of 2.5A. However, 
there are a few outlying points, both below (absorption) and 
above (emission) this band, which have | Wa | > 0.3A and 
which were selected for further investigation. In particular, 
there is a peak in the distribution of absorption features at 
6708A and emission peaks at 6717A and 6706A. 

The emission features at 6717A are due to [S 11] A6717, 
most probably originating in line-of-sight nebulosity in the 
LMC. This supposition is corroborated by the fact that the 
same stars also show strong Ha emission, as one would ex- 
pect. The occasional appearance of nebular emission lines in 
the 2dF data is not surprising because there is patchy nebu- 
losity across the LMC field, but the sky subtraction was done 
using the mean signal from a relatively small number of sky 
fibres chosen to lie in clear areas. However, the [S 11] A6717 
emission lines do help place limits on the acceptable values 
of the gaussian fwhm. Most of these have 2A < fwhm < 3A, 
and none has fwhm < 2A or fwhm > 4.5A. Consequently, all 
detections with fwhm < 2A were excluded from further con- 
sideration. In practice, no features rejected in this way were 
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close to the position of the lithium line. Since small spectral 
mismatches in the region around the lithium line can affect 
the wings of the detected lithium feature, as is the case with 
the strong line shown in Fig. 2c, an upper limit for the line 
width was set at fwhm < 5A. 

Another way to study these results is to consider 
Fig. 4 which shows histograms of the identifications plot- 
ted in Fig. 3 which lie in the range 2A < fwhm < 5A. The 
40 A wavelength range is divided into 2 A- wide bins. The 
three left-hand panels are for absorption features and the 
three right-hand panels are for emission features. The in- 
dividual histograms are for identifications with Wa > 0.4A, 
0.3 < Wa < 0.4A and Wa < 0.3A. The predominance of the 
peak at 6708A in the stronger absorption features is obvious 
(Fig. 4a): it contains eleven stars, one of which, KDM2832, 
has a very strong absorption line at 6707A and stands out 
clearly from the rest in Fig. 3. Its spectrum is plotted in 
Fig. 2c and shows a Li I A6707 that is strong enough to be 
seen even without the subtraction of the template spectrum. 
Of the remaining three features in Fig. 4a, one appears at 
669lA in star KDM 2379, another is the blend of two smaller 
features at 67I7A in star KDM 3389 and the third is the con- 
sequence of a bad pixel. Since none of the bins in Fig. 4a 
apart from the 6708A bin has more than one star, it is rea- 
sonable to assume that no more than one of the eleven stars 
in the 6708 A bin is there by chance. With only two out of 18 
bins having single stars, the likelihood that this is happening 
is only ~0.11. 

The same peak can be seen in the weaker lines of Fig. 4b, 
but here there are also significant numbers of lines in the two 
neighbouring bins at lower wavelengths. Nevertheless, the 
histogram certainly reflects a non-random distribution. It 
is unlikely that the detections shortward of 6708A are due 
to enhancements of individual metal lines because related 
metal lines are found in other wavelength bins as well (Barn- 
baum 1994) and are relatively weak in LMC stars. It is more 
likely that they are associated with a general spectrum char- 
acteristic which also produces the relatively high numbers of 
weak lines found (Fig. 4c) at the same two bins and also gives 
rise to a number of emission features at 6706A (Fig. 4e). This 
could be due to differences in the relative strengths of the 
12 C 12 C, 13 C 12 C, 12 C 14 N and 13 C 14 N molecular bands be- 
tween some stars and the mean spectrum. In contrast, very 
few very weak lines are found at 6708A. Nor are many emis- 
sion lines found at 6708 A. If, as is likely, the weak lines give 
a measure of how errors in the method generate lines, two 
of the ten lines in the 6708A bin in Fig. 4b are expected to 
have been caused in this way. If they do not, then the num- 
ber of false detections at 6708A (at 0.3 < Wa < 0.4A) could 
be as high as the number seen in the neighbouring bin at 
6706A, i.e., five out of the ten detections at 6708A could be 
false. 

Fig. 5 is an example of a J star with a modest detection 
of LilA6707 with W 6 707 = 0.48A. (This can be compared 
with the strongest lithium line detected in our sample, shown 
in Fig. 2c.) The superimposed lines show the quadratic 'con- 
tinuum' and the fitted gaussians. This particular example 
includes gaussians fitted to two broad shallow absorption 
features and two narrow emission features; all these have 
fwhm values that are inappropriate for real lines. Although 
very clear in this figure, the line at 6707 A cannot be seen 
directly in the original spectrum, but is manifested as a low- 
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Figure 4. Histograms of 'lines' detected in the residual spectra. 
Panels (a), (b) and (c) are for absorption 'lines' with Wa > 0.4A, 
0.3 < Wa < 0.4A and W A < 0.3A respectively; Panels (d), (e) and 
(f ) are for emission 'lines' in the same three ranges of Wa • 



ering of the maximum that usually arises between the groups 
of CN lines at 6702A and 6709A (see Fig. 1). Changing the 
'continuum' fit from quadratic to cubic or linear changes the 
set of selected lines slightly and generally reduces (cubic) or 
increases (linear) the equivalent widths. However, most lines 
selected as lithium are recovered in each case. 

At this point it is important to examine other possible 
sources of the features identified here. The CN bands coinci- 
dent with Li 1 A6707 are: Qi(28)+Q 2 (29) of the (12, 7) band, 
Q 2 (44) of the (6,2) band and Qi(22)+Pi 2 (13) of the (7,3) 
band (Davis & Phillips 1963). Within each of these CN sub- 
bands there are neighbouring bands which fall within the 
6690-6730A wavelength range and should be present at sim- 
ilar strengths. These are not seen in the residual spectra, so 
it is unlikely that the line detected at 6707A is due to CN. 
A similar argument applies to the Ri(17)+R,2(16)+R3(15) 
trio of lines belonging to the (1,4) Swan band of C2 (Phillips 
& Davis 1968), which appears at 6707A, because the neigh- 
bouring trios at 6703A and 6711 A are not seen. Other metal 
lines are unlikely to produce significant absorption near 
6707A (Barnbaum 1994). Another possible cause of spuri- 
ous spectral features is spectral mismatch due to a temper- 
ature difference between the star and the mean spectrum. A 
temperature difference of say 2000 K could produce a small 
apparent absorption at 6707A due mainly to altering the 
contribution from the (12, 7) band, but it would also ger- 
erate larger features at 6722A and 6728A. Since these are 
not seen, it is reasonable to assume that the 6707A fea- 
ture is not due to a temperature mismatch. The grouping of 
the stars in spectral subclasses helped to minimize any such 
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Figure 5. Normalized spectrum of an J-type carbon star 
(KDM2155) with a lithium line. The superimposed dotted and 
dashed lines are the sccond-oreder 'continuum' and the fitted 
gaussians respectively. 



problem. Finally, it is possible that in some cases the pres- 
ence of a circumstellar component in the 6707A absorption 
line can contaminate the photospheric feature. Although this 
possibility cannot be ruled out, it was not found to be a sig- 
nificant effect in Galactic carbon stars (Abia & Isern 2000) . 
Therefore, the observed absorption lines at 6707A seem to 
be correctly attributed to lithium. Ideally, however, high dis- 
persion spectra of these and some comparison stars would 
be needed to confirm the weakest identifications. 



3.3 Stars with detected Li I A6707 

Table la gives the final sample of carbon stars with absorp- 
tion features within 0.75A of Li I A6707.8, with W 6 707 >0.3A 
and with 2.0 A < fwhm < 5.0A. As seen in the previous sec- 
tion, features with W6707 > 0.4A are likely to be due to 
Li I A6707 at the 90 per cent confidence level, while those 
with 0.3 < W 6 707 < 0.4A are likely to be due to LilA6707 
at the 50-70 per cent confidence level. Column 1 gives 
the identification number of the star in the Kontizas et 
al. (2001) catalogue, column 2 gives the spectral group as- 
signed to each star (see Section 3.1), columns 3 to 6 give 
various photometric measurements which will be explained 
in Section 4, columns 7 and 8 give the W6707 (equivalent 
width) and fwhm measurements of the line at 6707A, col- 
umn 9 gives the j-index taken from Morgan et al. (2003) and 
column 10 provides cross references to Blanco, McCarthy 
& Blanco (1980): BMB, Blanco & McCarthy (1990): BM, 
Hughes (1989): SHV, and Westerlund et al. (1978): WORC. 

Table lb provides a list of LMC and SMC carbon stars 
with identification of LilA6707 recorded in the literature 
and makes a noteworthy comparison with Table la. SPLL95 
reported four Li-rich LMC carbon stars and a further 31 
carbon stars for which lithium could not be detected. Richer, 
Olander & Westerlund (1979) and Richer (1981) detected 
the lithium line in two out of 52 stars and two out of 25 
stars respectively. With one in each of the two pairs common 
to SPLL95, the total number of stars is six. SPLL95 also 
note two Li-rich SMC stars. These are both members of the 



Rebeirot, Azzopardi & Westerlund catalogue (Rebeirot et 
al. 1993) and both are J stars as seen in other 2dF data 
(Cannon et al., in preparation). They are also included in 
Table lb. 



3.4 Comparison with other work 

First there is the question of commonality between the 2dF- 
selected and previously known Li-rich stars. It turns out 
that not one of the six known Li-rich stars were included 
in the 2dF sample. Nor were any of the stars reported 
by SPLL95 as having no lithium detection included in the 
2dF sample. However, five of the other stars observed by 
Richer et al. (1979) (WORC 54 = KDM 1490, WORC 78 
ee KDM 1814, WORC 106 ee KDM 2518, WORC 192 ee 
KDM 5314, WORC 196 ee KDM 5426) and three of those 
observed by Richer (1981) (BMBBW31 ee KDM 2117, 
BMBBW41 ee KDM 2148, BMB BW 54 ee KDM 2197) are 
in the 2dF sample. Whereas LilA6707 was not noted by 
these authors as being present in any of these eight stars, it 
was found in one of the 2dF spectra (KDM 2197); none of 
the others has a line anywhere near the limit for inclusion 
in Table la. In fact, the line identified in KDM 2197 is at 
the lowest limits of both W6707 and fwhm for inclusion in 
the table. It is likely that this strength is below the limit for 
identification by Richer (1981). 

The equivalent widths, W6707, recorded in Table la 
range between 0.3A and 0.6A with one stronger line at 
1.4A. Torres- Peimbert and Wallerstein (1966) found simi- 
lar strengths for W6707 of ~0.4-0.5A in fourteen of the six- 
teen Galactic carbon stars in which they detected lithium. 
The other two were the well known super-Li-rich stars 
WZ Cas and WX Cyg. However, in the more recent stud- 
ies by SPLL95 in the Magellanic Clouds and by Boffin et al. 
(1993) in the Galaxy, significantly higher values W6707 have 
been reported. In fact, three of the four Li-rich LMC carbon 
stars described by SPLL95 are in the super-Li-rich category. 
Part of the discrepancy between levels of W6707 is due to the 
different methods used. SPLL95 used a pseudo continuum 
fitted near 6700A and did not attempt to correct their val- 
ues of W6707 for a CN contribution. It is worth stating that 
the We707 levels determined by Torres-Peimbert and Waller- 
stein (1966) and noted above as being similar to those of 
Table la were corrected for a contribution from CN, though 
through a method different to the one used here. Inspec- 
tion of Fig. 2c shows that this technique would give a larger 
equivalent width than the one quoted in Table la. Tests on 
the 2dF spectra of J stars imply differences of 0.4-0.5A be- 
tween the two methods. However, this level of discrepancy 
is not sufficient to explain why no super-Li-rich carbon stars 
were discovered in the 2dF sample. 



4 PHOTOMETRY 

Infrared photometry of sources in much of the LMC is 
now available through the Second Incremental Data Re- 
lease of the Two Micron All Sky Survey - (2-MASS) 
(Skrutskie et al. 1997) and is available on-line at 
http://irsa.ipac.caltech.edu/. 2-MASS JHK measurements 
were extracted for the stars in the 2dF sample by match- 
ing catalogue coordinates. Whenever 2-MASS data were 
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Table 1. Magellanic Cloud carbon stars with LilA6707 



Star 


Gp 


J 


K 


M bo i 




Ws707 


FWHM 


j-index Other 


a. 2dF stars 
















1310 


1 


11.37 


9.99 


-5.48 


3230 


0.49 


3.0 




1393 


4 


12.58 


10.74 


-4.46 


2670 


0.50 


4.0 


5.9 


1402 


1 


11.79 


9.86 


-5.32 


2580 


0.40 


3.0 


WORC 50 


1608 


4 


13.54 


12.46 


-3.34 


3740 


0.33 


2.5 


3.6 


1792 


1 


12.52 


10.97 


-4.37 


2990 


0.60 


4.0 




2155 


4 


12.61 


10.83 


-4.39 


2730 


0.48 


3.5 


6.0 


2197 


1 


11.58 


9.76 


-5.44 


2690 


0.30 


2.5 


BMB BW54 


2364 


4 


13.04 


11.20 


-4.00 


2670 


0.30 


3.5 


5.3 


2400 


1 


12.86 


11.34 


-4.02 


3030 


0.45 


4.0 




2406 


1 


11.37 


9.77 


-5.53 


2930 


0.35 


3.0 




2832 


4 


11.76 


10.34 


-5.09 


3170 


1.42 


5.0 


5.6 SHV 0515089-701859 


3344 


1 


11.24 


9.62 


-5.67 


2910 


0.31 


3.0 


WORC 138 


3495 


2 


12.53 


10.91 


-4.38 


2910 


0.44 


2.5 




3525 


1 


13.29 


12.10 


-3.56 


3530 


0.30 


3.0 




4404 


1 


12.53 


10.84 


-4.42 


2830 


0.35 


4.5 




4141 


4 


12.06 


10.56 


-4.81 


3060 


0.49 


3.5 


5.4 


4626 


4 


12.82 


10.84 


-4.33 


2530 


0.36 


5.0 


4.8 


4850 


2 


12.54 


11.19 


-4.31 


3270 


0.48 


3.0 


SHV 0534304-720850 


5455 


1 


12.65 


11.13 


-4.23 


3030 


0.59 


3.5 


BM 35-15 


b. Published stars 














1604 


4 


11.80 


10.36 


-5.00 


3140 


(4.43) 




WORC 65 


2041 


4 


13.30 


11.72 


-3.56 


2960 






BMB BW9 


2809 




12.33 


10.91 


-4.47 


3170 


(1.71) 




HV 5680 






11.98 


9.65 


-5.60 


2250 


(6.41) 




BMB R46 SHV 0521050-690415* 







11.48 


10.23 


-5.29 


3430 


(8.83) 




BMB BW89* 




2 


12.10 


9.35 


-6.22 


1980 






WORC 186* SHV 0535323-710211 




4 


12.64 


10.98 


-4.71 


2817 


(2.45) 




RAW 1296* 




4 


13.02 


11.29 


-4.37 


2740 


(1.55) 




RAW 1329* 



We707 and FWHM arc in A 

Wo707 values given in brackets were determined by a different method (see text) 

*: Gp represents an SC star (Richer & Frogel 1980) 

*: Gp 2 for WORC 186 is assumed from data by Richer ct al. (1979) 

*: BMB R46 is also F4488 as recorded by SPLL95 

*: RAW 1296 and RAW 1329 are SMC members 



unavailable, JK measurements were taken from the Deep 
Near-Infrared Survey of the southern sky - DENIS (Cioni 
ct al. 2000) which is available on-line from the CDS at Stras- 
bourg (http://cdsweb.u-strasbg.fr/denis.html). The magni- 
tudes from the DENIS database were corrected in the way 
described by Morgan et al. (2003) to bring them into align- 
ment with the 2-MASS data. 

The J and K magnitudes were subsequently used to 
derive bolometric magnitudes (M bo i) and effective temper- 
atures (T e ff) by applying the relations M bo i = K a — dm ± 
0.69 + 2.65(J - K) - 0.67(J - K)l and T eS = 7070/[(J - 
K) ± 0.88], derived by Wood, Bessell & Fox (1983) and 
Bessell, Wood & Lloyd Evans (1983), respectively. A dis- 
tance modulus value of 18.45 was adopted (Westerlund 
1997) and a mean reddening of E(J - K ) = 0.07 (Costa & 
Frogel 1996). Small deviations from the adopted value of 
the reddening do not significantly affect the derived values 
of T e g and M bo i- The resulting values of T ett - and JVf bo i are 
listed in Columns 5 and 6 of Table 1. It must be pointed out 
that the use of these parameterizations for Mboi and T e g 
may not be fully applicable to the specific stars studied here 
and must be treated with caution. Formal errors for Mboi 
and T c £f are ±0.2 mag and ±70 K. 



The bolometric magnitudes of the Li-rich stars identi- 
fied in the 2dF sample lie in the range —3.3 to —5.7 which 
includes fainter stars than the range of —4.6 to —5.7 occu- 
pied by the four Li-rich LMC stars identified by SPLL95. 
Similarly faint Li-rich carbon stars have been found in the 
SMC (SPLL95), in M31 (Brewer et al. 1996), and in the 
Galaxy where -3.5 < M bo i < -6.0 (Abia & Isern 1997). 

Figure 6 shows the location of the Li-rich carbon stars 
of Table la on the Hertzsprung-Russell Diagram. For com- 
parison the location of all carbon stars observed are shown 
on the same diagram. Their bolometric magnitudes and ef- 
fective temperatures were derived in the same way as those 
of the stars in Table 1, based on the data in Cannon et al. 
(2003: in preparation). The six Li-rich stars listed in Ta- 
ble lb are also included in this plot. The distribution of 
the Li-rich stars in Fig. 6 is clearly not the same as that 
of the normal carbon stars. Of the stars plotted as trian- 
gles, the two 'coolest' (WORC 186 and BMB R46) are both 
thought to be Mira variables and have the appropriate ex- 
treme colours, and the 'warmest' (BMB BW89) is actually 
an SC star (Richer & Frogel 1980) and perhaps should not be 
included. Two of the 2dF Li-rich stars are fainter and cooler 
than the others and lie near the tip of the red giant branch 
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Figure 6. A plot of M bo i against T e ff for LMC carbon stars. 
The circles are the Li-rich stars of Tabic la and the triangles 
are other Magellanic Cloud stars known to have strong Li I A6707 
(SPLL95). 
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Figure 7. Relative frequency of carbon stars in Mboi intervals 
of 0.5 mag. The left-hand panel is for N stars and the right-hand 
panel is for J stars. The open circles are for normal carbon stars 
and the solid squares are for the Li-rich stars in Table la. 



rather than on the AGB; the rest lie among the main body 
of AGB stars but seem to be more widely distributed and 
relatively uncommon at the peak of the general distribution. 
Fig. 7 shows the distribution of Li-rich stars with bolomet- 
ric magnitude (in 0.5-magnitude bins), separately for N and 
J types. For comparison, the distributions for normal car- 
bon stars are superimposed. The normal N and J stars stars 
have clear single peaks at M bo i = —4.75 and M bo i = —3.75 
respectively. Although the total number of J stars is small 
here, studies of a much larger sample show that the J stars 
are indeed fainter than the N stars (Morgan et al. 2003). 
The error bars shown in Fig. 7 are based on -Jn statistical 
errors. There would not seem to be a significant difference 



between the distributions of the bolometric magnitudes of 
the Li-rich and normal J stars. However, the Li-rich N stars 
are absent from the M bo i = —4.75 bin which is the peak of 
the normal N-star distribution. There are, in fact, no Li-rich 
N stars in the range —5.25 < Mboi < —4.5. 



5 DISCUSSION 
5.1 Statistics 

Of the 667 stars studied here, 53 were excluded as being 
of insufficient quality to allow the reliable non-detection (or 
detection) of LiiA6707 to a limit of W 67 07 ^ 0.3A. Of the 
remaining 614 stars, nineteen (i.e. 3.1 per cent) have de- 
tectable Li I A6707. Among the 614 stars there are 62 J stars, 
of which seven are Li-rich. This means that, to the detection 
limit of the present work, 11.3 per cent of the J stars are Li- 
rich compared with only 2.2 per cent of the non-J stars. So 
lithium enrichment seems to be about five times more fre- 
quent among J-type carbon stars. No super-Li-rich carbon 
stars were discovered in the sample of 614 stars. It is prob- 
ably true to say that none was discovered in the full sample 
of 667 stars because most of the extra 53 stars could have 
revealed a very strong lithium line with W6707 ~ 4A. 

These percentages are significantly lower than those 
recorded in the literature for the Galaxy and the Magel- 
lanic Clouds, although in the latter case in particular the 
samples are one or two orders of magnitude smaller than 
the present one. Boffin et al. (1993) found that in a sample 
of about 200 Galactic carbon stars ~12.5 per cent were Li- 
rich and 1.5 per cent were super-Li-rich. (The Boffin et al. 
(1993) results are updates of generally similar results based 
on smaller samples, presented in the earlier papers of Abia et 
al. (1991, 1993).) More recently, Abia & Isern (1997) made 
an extensive study of 44 Galactic carbon stars (28 of N type, 
11 of J type, and 5 SC stars) and detected 7 Li in 82 per cent 
of the J stars and in 15 per cent of the rest. In the Magellanic 
Clouds, SPLL95 found four Li-rich carbon stars out of a to- 
tal of 35 LMC carbon stars in their sample (11.4 per cent) 
and 2 out of 25 (8 per cent) for the SMC. 

Part of the difference in the frequencies of Li-rich carbon 
stars reported by the different authors may be related to 
actual differences between the carbon star populations in 
the parent galaxies (ages, masses, metallicities), but part of 
it may be due to the different approaches adopted for the 
lithium detection in the spectra, and, most importantly, the 
correction (or lack of it) for the CN contribution and the 
choice of continuum. The size of the sample considered in 
each case and its initial selection must also be taken into 
account. 



5.2 Spectral properties 

In order to gain some insight into the formation of the ob- 
served population of Li-enriched carbon stars, it is useful to 
investigate the possible existence of correlations between the 
equivalent widths W6707 and other properties of the stars, 
such as bolometric magnitude, effective temperature, and 
spectral group (Groups 1-4 as shown in Table la). Figs 8a 
and 8b show W6707 plotted against Mboi for the 2dF results 
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Figure 8. W6707 against Mboi- Panel (a) is for the stars of Ta- 
ble la (2dF results) and Panel (b) is for the stars of Table lb 
(SPLL95 data). The symbols are: J stars - solid squares; N stars 
- open squares; C2-weak stars - open circles; unknown types - 
open triangles. Note the different ordinate scales. 



in Table la and the data from SPLL95 (Table lb) respec- 
tively. The SC star, BMB 89, is considered to be too different 
in spectral type and is not shown. The two points in brackets 
are estimates of W6707 based on the Li-indices published by 
Richer et al. (1979) and Richer (1981). No statistically sig- 
nificant correlations can be found for the N stars in Fig. 8a, 
but the bimodality of the distribution can be seen again. Nor 
did Abia et al. (1993) find any correlation between bolomet- 
ric magnitude and lithium abundance for their sample of 
Galactic carbon stars. However, if only J-type carbon stars 
are taken into account, a different picture seems to emerge 
with W6707 appearing to increase with bolometric magni- 
tude. For the 2dF sample plotted in Fig. 8a, this conclusion 
is based on just one star and cannot be taken as more than 
tentative; but it is more striking for the SPLL95 sample 
plotted in Fig. 8b. As noted earlier, the apparent offset be- 
tween the new data and the SPLL95 results can be partly 
explained by the different continua adopted. In interpreting 
these diagrams it should be remembered that some of the 
stars, especially those in the SPLL95 sample, are variables 
and could change in Mboi and T c s; e.g., BMB R46 appears 
in Table lb with Mboi = —5.6 and is recorded at Mboi = —6.3 
by Whitelock & Feast (2000). With these considerations in 
mind, it is clear that a larger statistical sample is needed to 
confirm this result. 

So it would appear that J-type stars hold some special 
position in the Li-enriched carbon star sample: (i) it seems to 
be about 5 times more likely to get a Li-enriched carbon star 
if it is also of J-type and (ii) there seems to be a correlation 
between bolometric magnitude and equivalent width of the 
lithium line for these stars. It must also be emphasized that 
the highest equivalent width of the lithium line in the 2dF 
sample occurs in a J-type carbon star. 

W6707 does not seem to display any obvious correlation 
with either T e s or, for the J stars, with the j-index derived 
by Morgan et al. (2003) and listed in Table la. However, it 
must be remembered that the size of the sample of J-type 
Li-enriched stars is just seven. 



5.3 Theoretical models 

Given that for the first time the search for Li-enriched 
carbon stars has been conducted within such a large and 
homogeneous spectral sample, it is worth discussing the 
theoretical implications of the results presented here. Ven- 
tura et al. (1999) produced a self-consistent description of 
time-dependent mixing, overshooting and nuclear burning in 
AGB stars, and showed that there is a narrow range of stel- 
lar masses in which, despite the activation of the Cameron- 
Fowler mechanism, carbon can survive long enough in the 
convective envelope to keep the C/O ratio over unity as well 
as increasing 13 C/ 12 C, and that consequently these objects 
can be observed as Li-rich C stars for a non-negligible frac- 
tion of their AGB life. Their model can account for: 

(i) the very rare luminous Li-rich carbon star with Mboi = 
—6.5. No such star exists in either our or SPLL95's sample 
of LMC Li-rich carbon stars. 

(ii) Li-rich J-type carbon stars with Mboi — —5.75 to —5.95. 
The brightest J-type Li-rich star in our sample has Mboi = 
—5.1, which is more than half a magnitude fainter than re- 
quired, for its formation to be attributable to this mecha- 
nism. 

(iii) Li-rich non-J carbon stars for short periods during ther- 
mal pulses. Although this scenario could explain some of the 
N-type Li-rich carbon stars, it cannot explain all of them, as 
the mechanism can operate successfully only for short peri- 
ods. Not enough details are given by Ventura et al. (1999) to 
allow for a stronger statement. In any case, the faint J-type 
Li-rich carbon stars in our sample cannot be explained in 
this way. To summarize, it would seem that only very few, if 
any, of the Li-rich carbon stars of Table la can be accounted 
for by these models; the predicted overall distribution of Li- 
enriched stars is not a good match to the observations. It 
should be noted that the models which reproduce the charac- 
teristics of J-type stars suffer from another deficiency; they 
predict that these stars should show enhancements of the 
s-process elements because they are thermally pulsing AGB 
stars which suffer the third dredge-up, but no s-process en- 
hancements are seen in J-type carbon stars in the Galaxy. 
Thus, there is a strong possibility that J-type stars are not 
in the AGB phase and their chemical peculiarities are due 
to a different atmospheric processes. 

Travaglio et al. (2001) investigated the behaviour of 
AGB models - again in the framework of the HBB process 
- with masses in the range 4-6M and metal abundances 
between Z = 0.004 and Z = 0.02. They showed that the con- 
dition C/O > 1 together with high 7 Li (loge(Li) ~ 4) can 
only be achieved for the 4M Q Z=0.004 case, while much 
lower Li abundances (loge(Li) ~ 2) are found for the same 
mass, but for Z = 0.008. However, their models are far too 
bright (—7 < Mboi < —6), to be applicable for the cases of 
Table la. They do, however, mention that preliminary calcu- 
lations showed that lithium enrichment of AGB atmospheres 
through the HBB process can operate at lower masses (and 
give lower Mboi) for low metallicities (Z = 0.0001). This 
statement holds promise with regard to the faint Li-rich car- 
bon stars in the Magellanic Clouds, but one has to await for 
more detailed models. 

Abia & Isern (1997) found observational evidence that 
supports the existence of a deep mixing mechanism operat- 
ing in low-mass AGB stars that could be responsible for the 
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creation of Li-rich, J-type, low-mass (and low luminosity) 
carbon stars. This mechanism could account for the fainter 
J-type, Li-rich stars of Table la, but cannot explain the 
equally faint N-type, Li-rich carbon stars found. 



6 SUMMARY 

The results of the present study can be summarized as fol- 
lows: 

(1) The incidence of lithium enrichment among carbon stars 
in the LMC is much rarer than in the Galaxy, and about five 
times more frequent among J-type than among N-type car- 
bon stars: of the 614 stars studied here, 19 (3.1 per cent) 
have detectable Li I A6707. Among the 614 stars there are 
62 J-type stars, 7 (11.3 per cent) of which are among the 
Li-rich stars found here. No super-Li-rich carbon stars were 
found. 

(2) The equivalent widths of the LilA6707 line were mea- 
sured from spectra corrected in a statistical fashion for the 
contribution from (mostly) CN. The values derived are sim- 
ilar to the corrected equivalent widths for the lithium line in 
Galactic carbon stars derived by Torres-Peimbert & Waller- 
stein (1966) but significantly lower than the values derived 
by SPLL95 in the Magellanic Clouds (whose values were not 
corrected for CN contamination). 

(3) The bolometric magnitudes of the Li-rich carbon stars 
range between —3.3 and —5.7. The distribution in Mboi of 
the Li-rich N-type carbon stars seems to have two distinct 
peaks at ~ —4.25 and ~ —5.25 mag. 

(4) J stars hold some special position in the Li-enriched car- 
bon star sample, as they are about five times more likely 
to be Li-enriched. Also, there seems to be some correlation 
between bolometric magnitude and equivalent width of the 
lithium line for these stars. 

(5) Existing models of lithium enrichment via the hot bot- 
tom burning process fail to account for the observed prop- 
erties of the Li-enriched stars studied here. A different deep 
mixing mechanism operating in low-mass AGB stars seems 
to be required, especially for the faint N stars that do show 
Li-enrichment. 
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